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The  Partition  Wavenumber  in  Acoustic 
Backscattering  from  a  Two-Scale  Rough 
Surface  Described  by  a  Power-Law  Spectrum 

Jorge  C.  Novarini  and  Jerald  W.  Camthers 


Aboraa — Scatteriiig  tnm  the  ocean  bottom  Is  often  assumed 
to  be  controlled  by  two  spatial  scales:  the  iar|er  scale  associated 
with  reflections  from  plane  facets,  and  the  smaller  one  assoHated 
with  dWtase  scattering  ftom  height  variations.  Choosing  the 
wavenumber  for  this  partitioning  has  proven  to  be  important  bat 
troobiesome.  For  this  work,  scattering  data  are  simaiated  nsiag 
Hefanholtz-Kirchhoff  or  physicat  optica  theory  and  adectod  inpot 
geomorphology.  These  data  are  inverted  to  provide  rms  slope  of 
facets  and  rms  heights  of  small-scale  roaghnem  naing  a  afaapie 
two-scale  rooghness  model  introdneed  recently  (J.  W.  CamttMrs 
and  J.  C.  Novarini,  IEEE  J.  Octtmk  En^.,  voL  18,  pp.  108-188, 
1993).  Bottom  relief  is  described  by  power  spectra  of  the  powo* 
law  form,  and  the  bottom  is  assomed  to  be  impenetrable.  The 
work  introduces  a  new  criterion  for  effecting  this  partition  baaed 
on  setting  a  roughness  parameter  equal  to  unity.  The  crttcrioa  is 
shown  to  be  valid  for  the  cases  analyzed  baaed  on  the  ability  t/t 
die  inversion  model  to  recover  the  input  geomorphology. 


I.  Introduction 

The  two-scales  roughness  theory,  also  called  the 
composite-roughness  theory,  is  particularly  attractive 
because  it  separates  scattering  into  what  appears  to  be  its  two 
primary  mechanisms:  scattering  controlled  by  large  random 
facets,  and  scattering  controlled  by  smaller  scale  roughness. 
The  reported  advantage  of  this  iqjproach  is  that  its  tqjplicability 
can  span  a  wide  range  of  surface  roughness  since  it  allows 
for  the  selective  application  of  tqjpropriate  approximations 
in  scattering  theory.  While  this  is  certainly  true,  it  leads 
one  to  attempt  to  answer  the  troublesome  question  of  where 
to  set  the  scale  partition  based  on  the  requirement  for  the 
validity  of  the  qjproximations.  The  question  should  rather 
be  posed  in  terms  of:  “At  which  scales  is  the  scattering-like 
reflection  from  a  coUection  of  facets  and  at  which  scales  is 
the  scattering-like  diffraction  from  roughness?.” 

The  two-scale  roughness  (TSR)  iqiptoach  was  first  proposed 
by  Kur’yanov  [1]  and  further  developed  and  extended  by 
McDaniel  [2].  In  its  original  formulation  and  in  most  of 
its  ai^lications  [1H3],  TSR  scattering  models  resort  to  a 
combination  of  the  Helmholtz-Kirchhoff  or  physical  optics 
theory,  which  is  used  to  describe  scattering  from  the  large- 
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scale  surface,  with  a  Rayleigb-Rice  perturbatitHi  method  to 
account  for  scattering  at  the  small-scale  roughness.  McDaniel 
[4]  later  developed  a  TSR  scattering  model  based  entirely 
on  the  Kirchhoff  approximation,  and  derived  a  diffractive 
correction  to  relax  the  effects  of  tlK  selection  of  the  partition 
wavenumber  of  the  surface  spectrum.  Other  TSR  af^noaches, 
of  a  heuristic  nature,  that  allow  for  a  rapid  estimation  of 
the  scattering  strength  of  a  rough  surface  are  found  in  the 
literature.  Among  them  are  Brekovskikh  and  Lysanov  [S]  for 
scattering  from  the  sea  surface,  and  Martin  [6]  and  Ellis  and 
Growe  [7]  for  scattering  from  the  seafloor. 


n.  ApnjCABLE  Scattering  Models 

In  this  woric,  we  deal  with  a  TSR  model  originally  developed 
for  bistatic  bottom  scattering  strength  calculations  for  wide 
area  applications  [8],  the  bistatic  scattering  strength  model 
(BISSM).  An  upgrade  of  the  forward  scattering  calculation  of 
die  algorithm  has  been  tecendy  reported  [9].  The  term  BISSM 
has  been  ^plied  to  a  collection  of  programs  for  the  prediction 
of  the  scattering  strength  and  reverberation  from  the  seafloor. 
Here,  we  introduce  the  term  rough  facet  model  (RFM)  to  refer 
to  that  part  of  BISSM  that  is  the  mathematical  description  of 
the  scattering  mechanisms. 

The  application  of  the  rough  facet  model  to  sidescan  sonar 
images,  as  well  as  preliminary  applications  to  the  inverse 
problem  have  been  repented  elsewhm  [10],  [11].  The  attention 
of  this  work  will  be  focused  on  the  partition  wavenumber. 
The  ultimate  goal  is  to  allow  the  scattering  phenomenon 
itself  to  suggest  what  should  be  the  scale  partition.  RFM 
allows  inversion  of  the  scattering  problem,  and  thereby  the 
extraction  of  slt^  and  height  roughness  for  the  scattering 
surfaces.  Since  no  reliable  experimental  data  are  available 
(there  is  a  lack  of  scattering  measurements  accompanied 
by  a  detailed  description  of  the  bottom  roughness  at  the 
appropriate  scales),  scattering  data  are  simulated  using  a  high- 
fidelity,  2-D  scattering  model  and  selected  input  seafloor 
morphology.  This  high-fidelity  model,  which  we  refer  to  as 
the  benchmark  or  groundtrutii  model,  is  based  on  Mtmte  Carlo 
evaluation  of  the  Helmholtz-Kirchhoff  (H/K)  integral  over  an 
ensemble  of  surfaces  with  the  desired  statistical  properties. 
RFM  predictions  are  fitted  to  results  from  the  H/K  simulations, 
and  the  scale  partition  is  established  when  recovery  of  the 
ctHiect  geomoiphology  is  accomplished. 
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A.  The  Rough  Facet  Model 

The  RFM  algorithm  predicts  the  bistadc  scattering  strength 
of  a  directional  2-D  seafloor  for  arbitrary  source/receiver 
directions.  Here,  it  is  applied  to  monostatic  scattering  from 
2-D  isotropic  surfaces.  It  resorts  to  the  TSR  concept  in  the 
sense  that  the  problem  of  scattering  from  a  multiscale  rough 
surface  is  approximated  by  separating  the  scattering  problem 
into  fine-scale  and  microscale  roughness  (microroughness) 
effects.  Large-scale  roughness,  associated  with  long  surface 
wavelengths  present  in  standard  bathymetric  databases,  are 
handled  deterministically  and  are  excluded  from  this  analysis. 
That  is,  we  consider  here  only  the  roughness  within  the 
fooQ>rint  of  ensonification,  and  assume  that  any  tilt  of  the 
seafloor  under  this  footprint  can  be  removed  deterministically. 

The  scattering  coefficient  is  decomposed  into  two  com¬ 
ponents:  one  (component  mi)  mainly  controlled  by  diffuse 
scatter  in  the  incoherent  field,  and  the  other  (component 
m2)  mainly  controlled  by  a  facet  reflection  process  which 
for  individual  facets  is  partially  coherent,  but  collectively 
incoherent.  Accordingly,  the  total  bistatic  scattering  strength 
is  given  by 


mbt  =  mi  +  m2.  (1) 

The  adopted  methodology  is  somewhat  similar  to  the  treat¬ 
ment  followed  by  Brekovskikh  and  Lysanov  [5]  and  Martin 
[6]  (referred  to  as  “standard”  here),  in  which  the  mi  and 
m2  correspond  to  the  micro-  and  fine-scale  components, 
respectively.  However,  unlike  all  previous  work,  the  m2  term 
in  RFM  contains  information  on  microroughness  as  well  as 
fine-scale  slopes.  That  is,  the  major  distinction  between  FRM 
and  the  standard  model  is  the  distinction  between  rough  facets 
and  flat,  smooth  facets.  It  is  important  to  recognize  that  for  a 
random  distribution  of  flat,  smooth  facets  (a  broken  mirror), 
each  facet  reflects  coherently;  incoherence  is  introduced  by  the 
statistical  sum  over  the  distribution.  In  the  standard  treatment, 
as  well  as  in  the  methodology  adopted  here,  wide-angle 
scatter  is  produced  by  mi.  The  fine-scale  surface  gives  rise 
to  the  near-specularly  scattered  lobe  described  by  m2.  In 
contrast  to  the  standard  treatment,  however,  microroughness 
also  contributes  to  additional  loss  in  the  specular  lobe,  and 
provides  a  correction  to  account  for  diffractive  loss  in  RFM 
in  the  near-specular  direction. 

The  incoherent  component  is  based  on  Lambert’s  law.  It 
would  usually  include  an  empirically  determined  coefficient 
such  as  Mackenzie’s  coefficient  [12].  This  term,  however, 
is  ignored  here  because  we  limit  our  attention  to  the  near- 
specular  direction,  that  is,  in  the  monostatic  case,  near  normal. 
On  the  other  hand,  the  facet  scatter  part,  which  constitutes  the 
major  scattering  component  in  the  specular  lobe,  is  the  most 
sensitive  to  the  partition  wavenumber.  We,  therefore,  focus 
our  attention  on  m2. 

The  additional  coherence  or  scattering  loss  at  each  facet 
in  the  RFM  is  accounted  for  in  m2  with  a  multiplicative 
factor.  This  factor  was  chosen  to  be  the  so-called  “coherent- 
reflection  coefficient”  originally  derived  by  Eckart  in  the 
framework  of  an  H/K  formalism  [13].  It  accounts  for  the  low- 
flrequency  (i.e.,  small-scale)  scattering  by  the  microroughness 


in  the  local  specular  direction.  (By  “specular  direction,”  in 
general  we  mean  the  direction  a  ray  would  reflect  from  a 
hypothetical  flat  mean  plane  substituted  for  the  scattering 
patch.  The  additional  modifier  “local”  designates  the  specular 
direction  of  an  individual  fine-scale  facet.)  Therefore,  in  the 
RFM,  the  m2  component  of  the  scattering  coefficient  reduces 
to 


m2  =  M{k,  o^,S)-  F{6f,  6)  (2) 

where  M(  )  is  the  contribution  of  the  microroughness  and 
F(-)  is  the  contribution  of  the  fine-scale  slqpes.  The  additional 
parameters  are:  k,  the  acoustic  wavenumber;  o^,  the  standard 
deviation  of  the  microroughness;  £/,  the  standard  deviation  of 
the  fine-scale  slopes;  and  6,  the  angle  of  incidence  measured 
from  normal  (zenith  angle).  (We  also  define  the  slope  angle 
a/  =  tan“^6/.)  For  the  developments  of  Af(  )  and  F(-), 

Gaussian  height  and  Gaussian  slope  distributions  are  assumed, 
respectively. 

The  coherent  reflection  coefficient  is  given  by 

M  =  exp[-g^]  (3) 

where  g  is  the  well-known  roughness  parameter  [14]  restricted 
here  to  the  microroughness,  i.e., 

g^  =  4o^k^  cos^  0.  (4) 

The  function  F(-)  is  the  high-frequency  limit  of  the  H/K 
theory  as  derived  by  Eckart  [13],  Beckman  and  Spizzichino 
[14],  and  Brekhovskikh  and  Lysanov  [5].  (The  form  adopted 
here  is  identical  to  that  of  [S].)  For  the  present  q>pIication, 

F(-)  reduces  to 

where  Ro  is  the  Rayleigh  reflection  coefficient.  For  this  work, 
iio  =  1- 

In  [9],  we  introduced  RFM  in  a  heuristic  matmer.  Here,  in 
the  Appendix,  we  provide  a  more  formal  theoretical  foundation  ^ 

for  the  model. 

I 

B.  The  Benchmark  Model 

As  a  benchmark,  we  adc^t  here  a  model  which  consists  of 
numerically  evaluating  the  full  Helmholtz-Kirchhoff  integral  < 

with  a  minimum  number  of  tqqnoximations.  A  complete  | 

description  of  the  technique  can  be  found  in  [IS].  The  complex 
scattered  field  for  a  point  source  and  a  2-D  surface  is  calcu-  < 

lated.  No  Fraunhofer  or  Fresnel  iq)proximation  is  invoked,  and 
the  full  wavefront  curvature  is  kept  in  the  calculations.  The 
complex  pressure  at  the  point  r  is  given  by 

Ri  •  R2  \  Ri  R2 ) 

where  Ri=r'- fi^urce  and  R2  =  Rncemt-r',  and  r'  is  the 
position  vector  of  a  point  on  the  random  surface.  The  integra¬ 
tion  indicated  is  over  the  2-D  scattering  surface  S',  and  n  is 


tan*0 


26j 


(5) 


202 


IEEE  JOURNAL  OF  OCEANIC  ENGINEERING.  VOL.  19,  NO.  2.  APRIL  1994 


Fig.  t.  Geometry  of  the  simulated  scattering  experiments. 

the  normal  to  the  surface  element  dS*.  The  function  D{x,  y) 
represents  the  ensonification  function  over  the  surface.  Here, 
for  the  monostatic  case,  jiisourcel  =  IHreceiverl  =  R- 

There  are  two  advantages  for  using  this  model  as 
groundtruth:  the  H/K  formalism  has  been  shown  to  be  highly 
accurate  in  predicting  near-specular  scattering  ftom  multiscale 
surfaces  [16],  and  the  RFM  model  is  based  entirely  on  the 
H/K  approximation. 

III.  Numerical  Simulations 

Although  this  is  a  purely  numerical  study,  a  practical  sce¬ 
nario  based  on  a  high-frequency  multibeam,  swath  bathymetric 
sonar  system  is  adopted  here  because  it  will  allow  easy  testing 
of  the  results,  and  because  we  envision  that  system  as  a 
major  application  of  these  results  once  they  are  validated. 
However,  the  scenario  is  also  applicable  to  the  low-frequency 
active  sonar  system  used  in  the  Acoustic  Reverberation  Special 
Research  Program  (ARSRP)  of  the  Office  of  Naval  Research 
(ONR).  Data  taken  during  a  recent  experiment  of  the  ARSRP 
will  be  useful  in  future  validations  of  results  obtained  here. 
In  particular,  the  lowest  frequency  (250  Hz)  chosen  here  is 
applicable  specifically  to  that  data.  The  highest  frequency  (2 
kHz),  while  still  lower  than  the  range  of  standard  shipbome 
multibeam  sonars,  is  within  the  range  of  new  or  planned  towed 
swath  systems  that  incorporate  both  sidescan  imaging  and 
bathymetric  ciqiabilities. 

A.  Description  of  the  Numerical  Experiments 

Fig.  1  illustrates  the  geometry  of  the  simulated  experiment 
Shown  in  the  figure  is  the  representation  of  a  typical  transverse 
beam  of  the  m^tibeam  sonar.  The  water  depth  is  chosen  to  be 
2000  m  and  the  beamwidth  is  2°.  Zenith  angles  are  set  at  3-43° 
in  2°  increments.  Four  acoustic  frequencies  ate  investigated: 
250  Hz,  500  Hz,  1  kHz,  and  2  kHz.  A  constant  velocity  of 
15(X)  m/s  is  assumed.  The  use  of  an  isovelocity  sound  velocity 
profile  is  not  a  problem  for  the  near-vertical  angles  used  in 
the  study. 

For  each  ensonified  region,  the  pressure  was  calculated 
from  (6)  assuming  a  cosine-squared  ensonification  function  to 
minimize  sidelobes;  the  scattc^  intensity  is  then  obtained  as 


the  ensemble  average  of  pp*  over  the  ensemble  of  subregions. 
The  backscattering  coefficient  is  defined  as 

BSS  =  101og((pp’>flV^/i))  (7) 

where  A  is  the  area  of  the  ensonified  region  (weighted  by 
the  ensonification  pattern)  and  li  (=  pl/R^)  is  the  incident 
intensity.  We  normalize  the  source  pressure  (p«)  at  1  m  to  1 
/xPa. 

B.  Simulated  Two-Dimensional  Surfaces 

Ensembles  of  2-D  (i.e.,  surface  heights  are  functions  of 
X  and  y)  randomly  rough  surfaces  were  generated  using  a 
filtering  technique  developed  by  Caruthers  and  Novarini  [17], 
[18].  This  technique  consists  of  smoothing  a  2-D  array  of 
uncorrelated  random  numbers  with  a  filter  whose  transfer 
function  is  proportional  to  the  square  root  of  the  power 
spectrum  of  the  surface  elevations.  The  amplitude  spectrum 
that  is  most  often  used  to  describe  seafloor  roughness  is  of 
the  power-law  form 

A(A')  =  oA'-*’  (8) 

and  the  surfaces  are  assumed  to  be  isotropic.  K  is  the 
wavenumber  of  the  surface  roughness  (in  cycles/meters).  Note 
that  the  corresponding  power  spectrum  is  W{K)  = 

It  should  be  noted  ttud,  in  the  RFM  scattering  model,  the 
two  surface  parameters  (the  rms  height  of  the  microscale 
surface  and  the  rms  slope  of  the  fine-scale  surface)  are  band- 
limited  quantities.  The  variance  of  the  microroughness  for  a 
2-D  isotropic  surfaces  is  given  by 

al  =  271  W{K)K  dK.  (9) 

JKc 

The  variance  of  the  slopes  of  the  fine-scale  surface  is  given  by 

fK, 

6j  =  2x  W{K)K^dK.  (10) 

Jk, 

In  (9)  and  (10),  Ki  and  indicate  the  low-  and  high- 
frequency  wavenumbers  present  on  the  surface,  respectively. 
These  wavenumbers  are  ttetermined  by  the  size  of  the  fooqnint 
and  the  grid  spacing,  respectively.  Ke  denotes  the  partition 
(cutoff)  wavenumber.  When  the  partition  wavenumber  is  al¬ 
lowed  to  be  frequency  dependent,  then  both  and  ate 
functions  of  the  acoustic  frequency.  This  point  is  made  to  draw 
attention  to  the  fact  that  F(-)  of  (5)  is  frequency  dependent, 
contrary  to  the  usual  presumption. 

Table  I  shows  the  spectral  parameters  a  and  b  for  the 
surfaces  used  in  this  work.  Given  in  the  table  ate  ot,  the  total 
rms  height  (i.e.,  allowing  all  the  deformation  of  the  surface  to 
be  microroughness)  and  Sr,  the  total  rms  slqie  (i.e.,  allowing 
all  the  deformation  of  the  surface  to  be  slope  variance).  Also 
given  in  the  table  is  the  slqie  angle  ut-  llie  dimensions  for 
each  of  the  surfaces  is  256  x  256  m  (1024  x  1024  grid  points, 
with  a  grid  spacing  of  0.25  m).  The  grid  spacing  of  0.25  m 
bounded  the  analysis  in  frequency  and  an^e  to  2  kHz  and 
43°,  respectively,  owing  to  the  existence  of  grid  resonances 
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TABLE  I 

Parameters  of  the  Power-Law  Sfectra  for  the  Simulated  Surfaces 


Sutface 

Label 

Spectral  Paraameters 
a  (m^)  b 

nns  values 

ir  “T  (deg) 

OF 

2.32  xKT* 

1.75 

1.80 

0.594 

30.7 

01 

1.86x10-^ 

1,75 

1.45 

0.477 

25.5 

OD 

6.00X  10-^ 

2.27 

4.47 

0.464 

24.9 

OI 

4.80X  10-5 

2.27 

3.58 

0.372 

20.4 

OJ 

3.0X  10-5 

2.27 

2.24 

0.233 

13.1 

beyond  these  limits.  For  Kiichhoff  validity,  however,  we  will 
not  attempt  to  fit  data  beyond  35°. 

Each  surface  realization  was  split  into  a  number  of  en- 
sonified  regions  (sonar  footprints)  for  ensemble  averaging 
purposes.  The  number  of  subregions  depends  on  the  geometry 
of  the  simulation  and  the  zenith  angle  of  incidence,  and  varies 
between  9  to  2  as  the  angle  varies  between  3  and  43°. 

IV.  Results 

To  investigate  the  spectral  partition,  predictions  from  the 
RFM  model  were  fitted  to  the  simulated  data  from  the 
groundtruth  H/K  model.  These  fits  established  values  for 
(7 ft  and  6f  which  were  required  to  be  consistent  with  the 
simulated  surface  geomorphology  and  were  constrained  by 
(9)  and  (10).  Success  was  achieved  when  a  value  of  Kc  could 
be  established  that  satisfied  the  constraints.  As  mentioned 
earlier,  the  partition  wavenumber  is  usually  established  based 
on  either  the  formal  limits  of  the  approximations  invoked  [19] 
or  upon  a  partition  factor  that  is  the  ratio  of  the  partition 
wavelength  for  the  surface  Lc  to  the  acoustic  wavelength  la, 
that  is,  X  =  Le/la  [20].  The  cutoff  wavenumber  is,  therefore, 
Kc  =  k/{2‘Kx). 

Here,  we  hypothesize  that  the  scattering  phenomenon  itself 
will  establish  a  facet  size  such  that  the  microroughness  that 
exists  upon  it  leads  to  g,t  ss  1.  Larger  regions  contain  too 
much  roughness  to  produce  a  coherent  facet  reflection  and, 
therefore,  cannot  be  a  facet.  Smaller  regions  are  not  the  full 
extent  of  a  facet  that  can  cause  a  coherent  reflection.  Therefore, 
die  partition  wavenumber  can  be  determined  from  (4)  and  (9) 
and  6f  determined  from  (10). 

For  simplicity,  we  will  approximate  for  all  angles 
(0-43°)  by  its  value  at  0  =  0°,  i.e.,  =  4<r*A:*  »  1. 

This  could  lead  to  some  error  at  the  wider  angles.  However, 
we  do  not  press  the  issue  for  precision  because  other  sources 
of  error  exist  at  the  wider  angles.  The  sources  of  error  at  witte 
angles  will  be  discussed  in  the  next  several  paragr^s  as  we 
quantify  the  amount  of  error. 

Adopting  the  aforementioned  criterion  (i.e.,  to  assume  that 
=  I  Btiti  determine  the  resulting  Kc  and  £/),  we  proceed 
to  compare  the  scattered  field  firom  surface  OD  {b  =  2.27)  at 
S(X)  Hz,  with  predictions  from  the  RFM.  Fig.  2(a)  illustrates 
the  effect  in  (2)  and  (5)  of  varying  the  slope  of  the  spectrum 
(6)  in  computing  6f  from  (10).  Note,  first,  that  for  this  case, 
scattering  described  by  RFM  is  strongly  dependent  on  the 
spectrum  slope.  And  second,  the  best  agreement  occurs  at 
about  best  =  2.14.  This  departure  from  the  expected  b  = 
2.27  is  not  particularly  troublesome  because  spectral  analysis 


(a) 


0>) 

Hg.  2.  CranpaiisoDs  of  the  RFM  and  H/K  model  results  at  SOO  Hz  for 
surface  OD  for  different  values  of  spectral  slope  b  assuming  that  gfio-i-  (a) 
Backscattering  strength;  (b)  RMS  error. 

of  the  actual  surface  realizations  resulted  in  an  estimate 
of  spectral  slope  that  had  uncertainty  limits  of  2.10  and 
2.35.  For  practical  reasons,  the  surface  spectral  slopes  were 
estimated  from  the  overaU  surfaces  rather  than  averages  over 
the  ensonified  fooqnints  extracted  from  them. 

To  allow  for  a  more  detailed  examination  of  the  agreemenL 
we  establish  an  error  parameter.  To  that  end,  we  integrate 
the  difference  in  decibels  from  normal  incidence  out  to  some 
angle  6„  and  divide  by  6„  to  get  the  average  error  out  to  that 
angle.  In  digital  form,  this  becomes 

E  =  l/N*  SUM[abs(BSSRFM  -  BSSHK)](dB)  (11) 

which  provides  an  error  measure  from  normal  to  the  angular 
bin  N. 

We  examine  the  results  in  Fig.  2(a)  more  closely  by  applying 
(11)  and  display  these  details  in  Hg.  2(b).  Note  that  the  average 
error  is  increasing  with  angle  from  normal,  but  for  best  =  2.14, 
the  average  error  renuuns  below  0.5  dB.  It  might  be  noted  that 
the  point  errors  do  not  exceed  4  dB,  even  for  losses  of  33  dB. 
To  determine  if  we  can  fine-tune  the  value  of  g^to,  I^g-  3 
displays  the  enor  parameter  versus  angle  for  beat  =  2.14  and 
wiA  g^to  varying.  We  find  tiiaL  in  this  case,  =  1-0  iqrpears 
to  be  the  best  choice  to  two  significant  figures. 
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Hg.  3.  RMS  enor  between  RFM  and  H/K  model  leaula  for  suifoce  OD  for 
an  iteration  on  die  rooghnesa  parameter  of  the  micro-scale  surface  having 
established  the  spectral  slope  at  htes2.14  as  the  best  fit  for  gftt. 


TABLE  n 

BaNDUMTIED  RMS  RoUCMNESS  POR  SUMVtCE  OD  POR  gi/a«l  CORRESPONDINO 
TO  THE  Best  FtmNO  Parameters  a>i6.10'^  and  hew^^.U 


230  Hz 

500  Hz 

1  kHz 

2kHz 

<ru(m) 

0.47 

0.23 

0.12 

0.06 

ar  (deg) 

11.5 

13.1 

14.3 

15.3 

X 

4.9 

5.4 

5.9 

6.4 

Fig.  2(b)  and  3,  as  well  as  others  to  be  seen  later,  show 
eirors  beginning  to  appear  at  wider  angles.  Qeariy,  the  basic 
assumptions  of  near-specular  scattering  ate  breaking  down, 
and  we  introduced  an  obvious  error  by  using  instead  of 
g^,  but  there  are  odier  sources  of  error  introduced  by  isactical 
constraints  on  the  simulation.  These  include:  1)  at  tte  lower 
frequencies,  errors  are  introduced  by  difiraction  at  the  edges 
of  die  ensonified  region;  2)  at  the  higher  frequencies,  errors 
are  introduced  by  grid  resonances:  and  3)  statistical  errras 
introduce  the  limited  number  of  surface  realizations  (nine  near 
normal,  two  near  43° — although  not  as  bad  as  it  ^qiears  since 
there  are  many  correlation  lengths  in  the  surfaces  and  there 
are  more  at  wider  angles).  In  spite  of  diese  sources  of  error, 
the  results  are  very  good  out  to  35°  and  beyond  in  some  cases. 

To  test  the  robustness  of  the  critericm,  the  RFM  was  run  fm^ 
the  same  surface,  OD,  at  diree  odier  frequencies  (250  Hz,  1 
kHz,  and  2  kHz).  In  each  case,  diere  is  a  narrow  range  of  b 
values  that  agree  widi  the  H/K  curve  within  3  dB.  The  value 
6e,t  =  2.14  obtained  for  500  Hz  turned  out  to  be  a  good 
compromise  for  the  four  frequencies.  Hg.  4(a)-(c)  shows  the 
agreement  The  matching  between  the  RFM  prediction  and  the 
groundtruth  model  is  excellent  Thble  n  gives  the  mas  height  of 
the  microroughness  and  die  tins  slrqie  angle  of  the  fine-scale 
component  f<H'  each  frequency.  Also  shown  is  the  partition 
factor  X.  Notice  that  x  varies  significandy  and  consistendy 
over  the  three  octaves. 

The  same  procedure  of  kxddng  for  die  qiectral  paramo  b 
that  leads  to  the  best  agreemoit  for  =  1.0  was  repeated  for 
surface  OF  at  250  Hz.  (Note  that  surface  OF  has  a  significandy 
(Ufferent  total  variance  of  height  and  only  a  small  dififrience 
in  total  slope  variance  from  surface  OD.)  The  best  apeement 
(hart  =  1-75)  was  found  at  the  eiqiected  value  of  b  =  1.75. 
fr  Hg.  5,  die  effect  of  changing  is  shown  for  surface 


(>) 


(b) 


(c) 

Hg.  4.  Comparisons  of  the  RFM  and  H/K  model  results  for  surface  QD  for 
various  frequencies  having  adopted  fites2.14  and  (a)  2S0  Hz.;  (b)  1 

kHz;  and  (c)  2  kHz. 

OF.  As  luqipened  fOT  surface  OD,  the  ouqwt  of  RFM  is  quite 
sensitive  to  changes  in  p,,o-  Again,  the  best  agreement  occurs 
at  PmO  =  Ihble  in  lists  the  band-limited  roughness  ftn- 
surfaces  OD  and  OF,  at  250  Hz,  fm  different  p^o-  Ihe  partititMi 
facUns  (x)  are  also  shown.  It  can  be  seen  diat  x  ctHresptMiding 
to  the  be^  agreement  (p^o  =  1-0)  is  markedly  diffomit  in 
die  two  cases.  This  indicates  diat  die  numbn  of  acoustic 
waveloigdis  cannot  per  se  provide  a  criteritm  to  detmmiiie 
the  partition  wavenumber. 
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ARGU  OF  IHCIDDICl  (<l«g) 

Fig.  S.  Comparisons  of  RFM  and  H/K  model  results  for  surface  OF  for 
an  iteration  on  the  roughness  parameter  of  the  micro-scale  surface  having 
established  the  spectral  slope  at  t^=1.7S  as  the  best  lit  for  g^a=l. 


TABLE  m 

Bandumited  RMS  Microroughness  for  fs750  Hz  for  Numerical 
Realizations  of  Surfaces  of  (beti’°l.7S)  and  OD  (bM^2.l4) 


SURFACE  OF 

SURFACE  OD 

Xuo 

X 

<ru(m) 

oKdeg) 

X 

<»f(deg) 

0.2 

1.2 

0.21 

15.5 

2.4 

0.21 

13.3 

0.6 

2.5 

0.37 

12.1 

3.9 

0.37 

12.2 

1.0 

3.6 

0.47 

10.6 

4.9 

0.47 

IIJ 

1.4 

4.4 

0.56 

9.6 

5.7 

0.56 

11.0 

1.8 

5.3 

0.64 

8.9 

6.4 

0.64 

10.6 

To  jfuither  test  the  robustness  of  the  proposed  criterion 
for  the  partition  wavenumber  to  decreasing  roughness,  the 
model  was  ^lied  to  a  family  of  surfaces  (surface  01  and 
OJ)  that  were  generated  from  surface  OD  by  reducing  the 
surface  parameter  a  (i.e.,  decreasing  the  spectral  amplitude 
while  keeping  the  spectral  slope  6  constant).  Results  are  shown 
in  Fig.  6  for  a  frequency  of  250  Hz.  Table  IV  shows  the  rms 
values  of  the  fine-scale  slope  angle  and  the  rms  height  of 
the  microroughness  in  each  case.  In  all  cases,  the  criterion 
=  1.0  leads  to  the  best  agreement. 

As  a  final  application,  the  method  was  ai^lied  to  a  surface 
having  a  spectrum  of  the  type  describing  the  ARSRP  rectm- 
naissance  area  [18].  The  surface  is  labeled  01.  The  spectral 
parameters  are  listed  in  Table  I.  The  acoustic  frequency  is  250 
Hz,  and  the  partition  wavenumber  was  determined  according 
to  the  prtqiosed  criterion  9^0  =  1-  Fig-  7  shows  the  best 
agreement  with  the  RFM.  It  was  achieved  for  6e,t  =  1-85, 
i.e.,  5%  off  the  theoretical  value  of  1.75.  This  discrepancy 
is  well  within  the  error  in  the  spectral  slope  of  the  surface 
realizations.  The  resulting  microroughness  (<7^)  is  0.47  m  and 
the  rms  slope  angle  (a/)  is  12.3°. 

V.  CONCLUSKMfS 

The  partition  wavenumber  between  the  two  scales  of  rough¬ 
ness  in  conqxrsite-surface  scattering  the(»y  for  surfaces  having 
power-law  spectra,  typical  of  the  seailotMr,  has  been  inves¬ 
tigated  through  numerical  simulation.  Results  indicate  that  a 
criterion  based  (m  the  smoothness  of  the  microscale  component 


(») 


(b) 


Fig.  6.  Conqiaruoas  of  RFM  and  H/K  model  lesolix  for  (a)  nuface  01 
(bM=2.18)  and  (b)  sinftce  OJ  (ba>=2.26)  for  the  best  fit  at  gfio-t. 


TABLE  IV 

Best  Frr  Parameiers  and  Bandumiied  RMS  Mkrorouohness 
FOR  Surfaces  with  Different  Deoiees  of  Roughness  at  a 
Constant  Frequency  of2S0  Hz  (Curon'  Criterion:  gi/a’^l) 


luifaoe 

label 

ipecinl  paramelen 
a  (m*) 

band-ltd  nna  valoea 
(Tf,  Of  (deg) 

pMticioii 

X 

OF 

2.3x1(F 

1.75 

0.47 

lOJ 

3.6 

01 

L9x1(F 

1.85 

0.47 

12.3 

3.5 

OD 

6.0x10* 

2.14 

0.47 

IIJ 

4.9 

(H 

4.8x10-* 

Z18 

0.47 

10.6 

5.4 

OJ 

3.0x10-* 

2.26 

0.47 

8.9 

7.5 

may  be  adequate  to  establish  the  partition.  We  find  that  the 
partition  wavenumber  {qrpears  to  be  set  by  die  physics  of 
scattering  to  be  =  1.0  within  two  significant  figures. 

We  introduce  the  hypothesis  diat  a  facet  will  be  formed  such 
diat  its  roughness  is  ^iproaching  a  value  that  begins  to  destroy 
die  coherence  of  scatter  from  the  fac^  Larger  facets  ate  not 
allowed  because  coherent  reflection  is  not  allowed.  Smaller 
facets  are  not  allowed  because  thoe  still  could  be  a  cohemit 
reflecti(»  firmn  a  large  one.  Fbr  the  cases  analyzed  here,  the 
hypothesis  appears  to  be  valid. 

This  hypothesis  undoubtedly  has  limitadtms.  While  die 
partidon  factor  per  se  is  not  a  basis  for  the  establishment 
of  a  criterion,  it  will  likdy  set  a  small-scale  limit  to  slope 
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ANGLE  OF  INCIDENCE  (deg) 

Fig.  7.  Comparison  of  RFM  and  H/K  model  results  for  surface  01  for  the 
best  lit  at  6^1=1.85  and  gpa=l. 

scanering;  that  is,  the  cutoff  wavelength  cannot  be  an  atbitrar- 
iiy  small  fraction  of  an  acoustic  wavelength,  while  allowing 
9^0  be  near  unity.  Clearly,  the  scale  would  be  t(X)  small 
for  physical  optics  to  apply.  That  is,  we  cannot  expect  the 
approximations  leading  to  H/K  to  be  violated  too  seriously. 
(We  are  less  concerned  with  minor  violations.)  Or,  if  the 
surface  is  so  smooth  that  <  1  over  the  full  f(X)tprint  or 
many  Fresnel  zones,  then  we  would  likely  find  only  coherently 
reflecting  facets. 

There  are  two  implications  of  these  results  that  are  worth 
further  note:  coherence  loss  (i.e.,  exp[-</^o])  due  to  nti* 
croroughness  is  constant  and  equal  to  1/e  (-4.3  dB),  and 
the  previously  assumed  frequency  independence  of  F(  ),  as 
derived  in  [5],  [13],  [14],  and  others,  does  not  hold  for 
power-law  spectra.  The  implicit  frequency  dependence  of  F(-) 
through  the  slope  of  the  spectrum  of  the  surface  and  the  cutoff 
wavenumber  (Kc)  is  strong. 

Appendix 

Derivation  of  the  Rough  Facet  Model 

The  F(-)  factor  in  (2)  represents  scattering  due  to  the  fine- 
scale  slopes  (S/).  It  is  identical  to  (9.9.14)  in  Brekhovskikh 
and  Lysanov  (B/L)  [5]  for  the  high-Bequency  limit,  where 
they  use  the  symbol  V  to  represent  the  Rayleigh  reflection 
coefficient  (here  Rq)  and  are  only  concerned  with  the  large- 
scale  slopes  (^)  representing  regions  of  stationary  phase.  This 
V  was  taken  outside  the  integral  in  B/L  (9.7.4),  with  the 
assumption  that  the  Rayleigh  reflection  coefficient  changed 
slowly  with  angle  over  the  integral. 

We  are  also  concerned  with  roughness  on  the  facets,  and 
treat  each  facet  as  a  scatterer  of  the  form  given  in  B/L  (9.8.7). 
Following  the  same  procedure,  we  replace  with  from 
B/L  (9.8.7),  and  we  get  (in  our  notation) 

-.exp(-ff^)  fig. 

When  this  is  put  into  B/L  (9.9.14),  (2HS).  are  obtained. 
We  have  made  two  additional  assumptions:  the  height’s  dis¬ 
tribution  is  Gaussian,  and  exp(-p^)  is  also  slowly  varying 
with  angle  over  the  integral  in  B/L  (9.7.4).  In  effect,  all  we 
are  doing  is  replacing  the  Rayleigh  reflection  coefficient  for 


IEEE  JOURNAL  OF  OCEANIC  ENGINEERING,  VOL.  19.  NO.  2.  APRIL  1994 

each  specular  reflection  from  a  facet  with  the  Eckait  coherent 
reflection  coefficient  times  the  Rayleigh  reflection  coefficient. 

As  already  pointed  out  in  a  previous  woric  [9],  where  the 
rough  facet  model  was  introduced,  the  factor  exp  (-p^)  can  be 
identified  with  one  of  the  diffractive  corrections  in  McDaniel’s 
[4]  approach  to  scattering  from  composite  surfaces  using  the 
Kirchhoff  approximation  for  both  the  large  and  small  scales. 
In  that  development,  the  additional  approximation  mentioned 
above  is  also  required. 
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